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O . Abstract 

Q^' Standard model factorization can account for the large rate for 

\ rj'K^. The large rate for B — > rj'Xg is more problematic and requires a 

^ I "cocktail" solution. Individual contributions from factorization, b 

I and "intrinsic charm" fall significantly short of the observed rate. In contrast, 

the observed t]' yields are easily explained in models with Bijj sg) ~ 15%. 
Factorization can account for the entire semi-inclusive ry' yield without vio- 
lating rare decay constraints. Implications for CP violation are discussed. 
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I. INTRODUCTION 



The CLEO collaboration |I| has reported sizable production of fast rj' mesons, 

B{B r]'Xs) = (7.5 ± 1.5 ± 1.1) x 10"^ (2.0 < p^,' < 2.7 GeV) (1) 

and a correspondingly large exclusive rate, 

B{B^ 7]'K^) = (7.8+1^ ± 1.0) X 10-^ (2) 

There is no evidence for rj'K* modes in the inclusive analysis, with most events lying at 
large recoil mass, mx,- The experimental cut on prj' corresponds to mx^ < 2.35 GeV. There 
is also an exclusive branching fraction upper limit for rj production, 

B{B r]K^) < .8 X 10-^(90% c.L). (3) 

In the standard model (SM) factorization estimates for semi-inclusive and exclu- 

sive |^,|^ rare B decay branching ratios (i?r's) are typically of order 10~^ and 10~^, respec- 
tively, so it is tempting to speculate that the large t]' yields are due to the intervention of new 
flavor physics. A natural candidate to consider is enhanced b ^ sg chromomagnetic dipole 
operators which can explain several apparent discrepancies between SM expectations and 
experimental results in B decays: the low semileptonic Br and charm deficit [^p|,p!0|,^,p 



the kaon deficit ||T0[, and the low r(A;,)/r(i?(i) lifetime ratio |Tj|. All can be accounted for 
if B{b — * sg) ~ 10% — 15%, which is an interesting range because it is naturally associated 
with TeV scale dynamics for quark mass generation 0. 

In what follows we critically examine possible mechanisms for explaining the large rj' 
yields in the standard model and in models with enhanced b sg. We find that standard 
model factorization can, in principle, account for the exclusive rj' yield. For a liberal range 
of factorization model parameters B{B^ rj'K^) lies in the range (1.1 — 5.8) x 10"^ for 
rris = .2 GeV and (2.3 — 12.1) x 10^^ for rris = .1 GeV. It is worth noting in this regard 
that the recent values of rris obtained in lattice and QCD sum rule studies are quite low. 



e.g., m7(2 GeV) ^ 128 ± 18 MeV [y] and m7(2 GeV) ^ 100 ± 21 ± 10 MeV The 



inclusive yield is more problematic. For the same range of parameters SM factorization^ 
gives v'-^s) ~ (0.5 — 2.5) x 10~^ including the experimental cut, with the largest 

yields corresponding to a fairly limited region of parameter space. Of course there could be 
additional non-factorizable corrections, e.g., at the 10% — 30% level. 

Atwood and Soni (AS) |]TB[ have suggested that the large inclusive rate is connected 
to the standard model QCD penguins via the gluon anomaly, leading to the subprocess 
b — s> sg* f]'sg. Taking a constant ggrj' vertex form factor i7(0,0,m^,), extracted from 
J/\l' ^ 7]''-), they obtain agreement with Eq. (|l|). More recently, Hou and Tseng (HT) [l^ 
have argued that the factor as implicit in H should be running, which would lower AS's 
result by roughly a factor of 3. However, both AS and HT have overlooked the leading 
"^^'/(^^ ~ ^'^') dependence of the form factor, where q is the virtual gluon's momentum. 
Including this dependence nominally reduces AS's result to B{B v'-^s) ~ 1.6 x 10^^ 
including the experimental cut, more than an order of magnitude below what is observed. 

AS have estimated the contribution from decay of intermediate charmonia and from rj'—rjc 
mixing to be B{B —>■ rj'Xs) ~ 1.1 x 10~^, including the cut. Halperin and Zhitnitsky have 
suggested that the elevated glue content of the t]' gives it a large intrinsic charm component 
which simultaneously accounts for the exclusive and inclusive r]' yields via b ccs Vj'X. 
However, as already noted in |l^, their prediction of BiB v(K) ~ 2B{B rj'K*) is 
inconsistent with the absence of rj'K* modes in the inclusive analysis. This alone implies 
that the intrinsic charm r]' yields can not account for more than about a tenth of Eqs. 
(0) and (@).i Finally, as a by-product of their SM analysis of 6 ^ sgg Wyler and Simma 



^Throughout we obtain results for inclusive r/' production for charged B decays. We will update 
our results in the near future to include the average over charged and neutral B decays. The 
branching ratio results are not significantly modified and our conclusions remain the same. 

^It was recently suggested [|l^] that the anomalously large Br in Eq. (||) can be explained by 
color-octet cc contributions to r/' production, in analogy with the color-octet mechanism in the 



3 



have obtained an O{10~^) Br estimate for B ^ K + gluehall |jT9|], indicating that the Br 
for h sgg srj' should not exceed 0(10"^). In our view, although each of the above 
mechanisms falls short by itself, a standard model 'cocktail solution' for the inclusive rj' yield 
is still possible given the large hadronic uncertainties involved, with factorization providing 
the largest contribution. 

In contrast, we find that models with enhanced b ^ sg readily account for the observed 
inclusive yield in the factorization model alone, obtaining B{B v'-^s) ~ (1-0 — 9.0) x 
10~^. The actual value largely depends on the choice of weak phases in the dipole operator 
coefficients, and B ^ rj' form factors. Moreover, we have checked that the large rj' yields 
are compatible with the potentially restrictive CLEO upper limits on 5 — (pXg [^] and 
B (j)K for a significant range of parameters in the factorization approximation. 
The recoil spectrum increases with nix^-, in qualitative agreement with observation. In the 
case of the anomaly mediated process h Tj'sg, HT argued that the entire inclusive yield 
can be reproduced if b —>■ sg is enhanced. Taking into account the previously overlooked 
"^^'/(^^"'^v) dependence of the ggrj' form factor, but ignoring the subleading dependence 
of as in HT, we find B{B ri'X^) ~ (0.3 - 1.6) x 10"^ for B{b sg) ^ 15%. Although 
significantly smaller than the rates obtained in the factorization approach it is clear that 
the sum of the two can easily reproduce Eq. (|l]) if 6 — sg is enhanced. 

An important feature of quark model B decays pointed out by Lipkin is the inter- 
ference between spectator (6 suu) and non-spectator (6 sss) production of t]' and rj 
mesons. It is constructive for exclusive t]'K modes and destructive for exclusive rjK modes, 
leading to order of magnitude smaller rates for the latter, in accord with the recent measure- 
ments. We have checked this explicitly in the factorization approximation for the standard 
model and for enhanced b sg. The constructive interference can also account for the larger 



NRQCD description of charmonium production. However, the cc pair must be far off shell inside 
of the r/', invahdating the use of the NRQCD formalism. 
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rates for rj'K versus i^vr, (pK. Lipkin has noted that the same interference pattern persists 
generally for orbital parity even final states, e.g., constructive for rj'K** and destructive for 
rjK**, while the reverse is true for orbital parity odd final states, e.g., destructive for rj'K* 
and constructive for rjK*. This most likely can explain the absence of rj'K*^ events in the 
inclusive analysis. We leave a factorization model analysis of B ^ rj'K* for future work. 

The weak phases associated with new physics contributions to the dipole operator co- 
efficients generally lead to large CP violating asymmetries in charmless B decays in the 
factorization approximation [^] when strong phases in the penguin amplitudes from cc 



rescattering [0,^ are taken into account. For SM factorization we obtain CP asymmetries 



in the range 4.5% — 8% for 5^ rj'K^ and 2% — 5% for B fj'Xg with, as usual, the larger 
asymmetries corresponding to smaller Br's. Including enhanced b ^ sg can give asymme- 
tries of order 10% (exclusive) and 5% (inclusive) for Br's in the observed range. These 
are not dramatically larger than the SM asymmetries, unlike the 'pure penguin' processes 
B"^ —>■ (f)K^ and B —>■ K^n^ for which large asymmetries are possible compared to < 1% in 
the standard model. We note that enhanced b ^ sg can also lead to large isospin violating 
asymmetries in radiative B decays [Q. 

In the next section we discuss t]' production in the standard factorization approach. 
In Section 3 we focus on the gluon anomaly subprocess and comment on other 'exotic' 
mechanisms. We conclude with discussion of our results in Section 4. 

II. T]' AND T] PRODUCTION IN THE FACTORIZATION APPROACH. 

The effective Hamiltonian for non-leptonic charmless b ^ s transitions is given belowi 



^The sign of the operators On, O'n corresponds to use of —igsT^ for the quark - gluon vertex 
Feynman rule , leading to destructive interference in the standard model between penguin and 
chromomagnetic dipole operator contributions to charmless b sqq decays. 
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n 



4G, 2 11 



lly 

i=l i=3 



03(5) = s'y^.LbJ^lYHR)'! , 04(6) = s„7^L6^^g/37^L(i?)g„ , 

3 3 
^7(9) = 2^7M^&IZe'?^7''^(^)9 , 08(10) = -Sa7M^^/3 II e5g;37''L(i?)g, 



Ou = ^m,{fi)sa,,Rt%G^'' , 0[, = ^mt{fi)sa,,Rt%G>^'' , (4) 

where the standard notations are used and q is summed over u, d, s. 

In our numerical analysis we take as{Mz) = .117, rrit = 174 GeV, rrif, = 5.0 GeV, 
irtc = 1.63 GeV (HQET parameter Ai = —.2), and = .1 to .2 GeV. For KM matrix 
entries we take VtsV^*^ = .037, and Wolfenstein parameters p = —.15, rj = .33 which are in 
the range favored by the SM analysis of Ali and London [^. We use the Next-to-Leading 
Order (NLO) scheme-independent effective Wilson coefficients p9|J^,|^ at /x = mf,, obtaining 

ci = -0.306, C2 = 1.15, C3 = 0.017 - 0.37Ps, C4 = -0.037 + l.llP^, 
C5 = 0.01 - 0.37Ps, C6 = -0.045 + l.llP^, C7 = 1.03 x 10"^ + 0.33Pe, 

Cg = 3.82 X 10"^ cg = -0.01 + 0.33Pe, Cw = 1.97 x 10^1 (5) 

The P,,, are given by P, = ^(f - = _ GK, /i, g^)), where 

G(m, /i, g2) = -4 /" x(l - x) In ~ ""^^^ ~ ""^^' (ia;. (6) 

^ 

g is the momentum of the virtual gluon in the penguin diagram. G becomes imaginary when 



q > Ami ^ strong phase is generated P4| , ^ 



The Leading Order (LO) dipole coefficient is c^™(m5) —.295, which corresponds to 
B{h sg) ~ .2%. For enhanced dipoles we'll consider E{h —>■ sg) ~ 15%, corresponding 
to@ (|ciip + |c']^iP)^ ~ 2.26. In obtaining Pr's we consider two total B decay widths: the 



^This value is determined using the total NLO standard model, inclusive b decay width, see next 
footnote, and the LO expression for B{b — > sg). 
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NLO SM value! for the above inputs, Tsm ^ 3.3 x IQ-^^ Q^y^ ^nd 3.88 x IQ-^^ Q^y 
B{b — i> sg) ~ 15%. The dipole coefficients are parametrized as 



Cll 



\Cu\e 



^11 



-\Cn\e 



(7) 



Generally, large direct CP violation effects result from interference between new amplitudes 
containing the weak phases 6*11, 9[-^ and standard penguin amplitudes containing strong 
phases from the G function. 

In what follows the dipole operators are included by allowing the off-shell gluon to turn 
into a quark-antiquark pair. For On this leads to the effective Hamiltonian 

.Gf Ci 



Till 



E 

q=u,d,s 



v/2 vrg' 



:VtbV*Gnmb{lJ,)s{ps)ai,^Rt%{pb)q{p2h'''taq{pi)q'', 



(8) 



where = [p^ — PsY ■ In the standard model Hn interferes destructively with the penguin 
amplitudes. In practice, we take = gs{iTLb)gs{q)/4:TT, identifying one factor of gs with the 
virtuality of the gluon. This procedure models the fact that in approaching the 'on-shell' 
gluon limit Tin becomes increasingly non-local so that gluon fragmentation must take place 
at non-perturbative size scales. Performing Fiertz rearrangement yields! 

2N, 



n 
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+ E 



q=u,d,s 



47rg2 



[2mbSa1^,Lqpqo,'YLhfii - ArribSaRqpqa'Lhpi + 2msSa-ff,Rqi3qa'l^R^i3' - 
ArUsSaLqpqa'Rbp' + {pb + Ps)f, s^'-f^Lq^q^iRbp' + SaRqpqa'l^Rhp' + 

iSaa'^'RqpqanuRbp' - iSaluLq^q^'a^" Rhpi | 

The effective Hamiltonian Ti'^^ corresponding to O'n is obtained by substituting C 



(9) 
C'n 



and L i? in the above. In all of our factorization model results we take 



^This includes the NLO corrected b ^ c and b ^ u transitions [^0|, taking into account 0{as(^) 
corrections to the LO penguin contributions to 6 — > ccs, and 0{l/ml) HQET corrections The 
charmless b ^ s transitions and On have not been included but their total correction is ~ 1%. 

^This result differs from the result in ||3|. 
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A. Semi-inclusive rj ^ rj production. 

Factorization model contributions to il^'^'Xs involve two classes of amplitudes 

distinguished by their hadronization pattern, which we loosely refer to as "two-body" and 
"three-body" decays.0 The matrix element is of the form 

M = M^^ + M^^ (10) 

so that 

|M|2 = IM^^I^ + |M^^|2 + M^''^M^^ + M^''^M^\ (11) 

In the two-body decays an t]' or rj is formed from an ss pair via the subprocess b sss, 
which in the parton model looks like b r]^'^s. Applying the Gordon identity and Dirac 
equation the matrix element for rj^'^ production can be written in the simplified form 

M^' = (f {Xs\sLb\B) + (f {Xs\sRb\B). (12) 

For 1]' production the coefficients depend on the matrix elements 

{V'\si^l5s\0) = -2p^,^(Focos^ - V2F8sin^) (13a) 

{V'\s75s\0) = t^^Fssine (13b) 

(V|m7^75u|0) = -zp^,^(Focos^ + -^Fgsin^), (13c) 

where Fq and Fg are the decay constants for the SU{3)p singlet and octet axial vector 
currents, and the anomaly has been taken into account in (?7'|s75s|0) following For the 
7] the corresponding matrix elements are 

(r/|s7^75s|0) = zp^^(Fosine + v^Fgcos^) (14a) 



similar treatment of — > KX is presented in ||5i 
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{r^\s-f,s\0) = ^Fgcos^ (14b) 

{r]\uYj,u\0) = zp^^{Fosm9~ l^FsCosO). (14c) 

Phenomenological fits in P2| give Fq ~ Fg ~ F^, where = 132/^2 MeV, and an f]' — Tj 
mixing angle 6'^^' ~ —17°. Chiral perturbation theory favors Fg ^ 1.25F^ and On ~ —21" 
33| , |32[| . fir's are given for both choices. We make the usual quark model approximation, 



taking equal momenta for the two s quarks in the r^*^'-*, which gives 

= ml/2- ml,, /4 + ml/2. (15) 
For simplicity we take the b quark and spectator at rest so that the rj^''^ energy is fixed. 



Taking their momenta into account as in Ref. would not have a significant impact on 
the overall Br but would smear the "two-body" rj^'^ spectrum by a few hundred MeV. 

SM b Tj's Br's can be found in Table |. Associated CP asymmetries are ~ +1%. 
The large dependence is traceable to the l/m<j factor in {r]'\s'j5s\0) which, as we'll see 
shortly, can lead to large exclusive fir's in the SM. The corresponding b —>■ rjs fir's range 
from (3.4 — 9.0) x 10^^, with CP asymmetries of ~ —1%. Enhanced b sg can increase 
these contributions by a factor of 2 to 3, but the dominant contribution will come from the 
"three-body" decays. 

The "three-body" decays involve hadronization of the spectator into the i]^'^ via the 
subprocess b — *■ suu. In the parton model this is interpreted as fi^ rj^'^su, with px = 
Ps +Pu- Using the equations of motion and again, for simplicity, taking the b quark at resti, 
the matrix element can be simplified to the form 

M'" = Ci'^^iXsuls^^LulB) + if {Xsu\sRu\B) 

Hf\x,^\s^^Ru\B) + if {X,^\sLu\B) (16) 



^There are in general two additional terms which vanish in this limit, being proportional to 
p^p^sa^jjRu and p^p^sa^^Lu following application of the equations of motion(EOM). 



The matrix elements entering the are parameterized as ||35||: 



(r/W|^7'^6|i?-) = {ps+p,,,rFf\mlJ + {ps - P,^,Y P'\m\^) (17a) 
{ri^'^\uh\B-) = -ml,,)Ff{mj,J+mlj^'\mlJ (17b) 

(V \ucr^ b\B ) = -I —{PBPnC) -PsPjw)' (l^c) 



where, 



;,'"(p.) = '"''"^;'"^- |Ff(/) - (18) 

With a simple monopole parameterization assumed for the form factors they are determined 
by fixing Fi~ and Fq+ at zero momentum transfer, where they must be equal. It appears, 
based on various models in the literature p5|,r7|,|5B|,p7|] , that a reasonably broad range to 



consider for all B to light pseudoscalar form factors is Fi-(O) = .25 — .5. 

To gain insight into the factorization approach it is useful to study the quark level NLO 
Dalitz plots for b sqq shown in Fig. (P, where = {pb — VsY and rn? = {ps + Pq)"^. 
For simplicity we only consider h —>■ sdd, which does not contain tree-level contributions. 
The same conclusions can be drawn from the Dalitz plots for b suu and b sss. We 
have blown up the regions relevant to factorization for two-body or semi-inclusive "quasi" 
two-body decays, i.e., ~ and m < 2 GeV . It is clear why models with enhanced 

b —>■ sg readily evade constraints from rare B decays in the factorization model. Although the 
charmless Br is an order of magnitude larger than in the SM it is peaked at low g^, i.e., the 
"on-shell" gluon limit.i In the region relevant to rare decays, e.g., B Ktc, Kcp, (pXg, etc.. 



^The singular low q"^ dependence of 6 — > sqq originates from the pure dipole contributions, e.g., 
0{cii). Although these contributions come in at 0{a'^) they are clearly important in the case 
of enhanced dipole coefficients and should be included, particularly since we will be interested in 
comparison of inclusive quark level and factorization model Dalitz plots. This dependence will be 
canceled in the total charmless Br by 0{as) corrections to 6 ^ sg. 
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the pure dipole amplitude is only of same order as in the standard model. This means that in 
general large interference is possible between the penguin and dipole amplitudes. Generally 
this will lead to large CP violation ||23| because at ~ mf/2 the penguin amplitude should 



develop large strong phases p^ , |25| , and the dipole amplitude can have arbitrarily large weak 
phases.0 

Three-body rj'su Dalitz plots are shown in Fig. @). Note that the 'quark level' 

definition for q is used, i.e., pb — Ps rather than ps — Ps- The plots do loosely resemble the 
corresponding inclusive plots in Fig. (|l]). The regions surviving the cut < 2.35 GeV 
are blown up. The corresponding recoil spectra, shown in Fig. @, show a characteristic rise 
with mx^, like the measured spectrum Including the two-body spectra and interference 
between two-body and three-body rj' production (2b-3b interference), particularly the de- 
structive interference expected for exclusive rj'K* production |2^, should improve agreement 
with the shape of the observed spectrum at low niXs- 

To address the question of whether 2b-3b interference can have a significant impact on 
the semi-inclusive r]^'^ yields we have blown up the regions of Figs. (Jlja) and (|T]c) which are 
relevant for three-body t]' production, with experimental cut, in Fig. (^. The three-body 
amplitude can interfere with two-body r]^'^ production when q ~ and mx^ ~ 1 GeV. 
A glance at Fig. (^) will convince the reader that in the SM the relative contribution of 
this region to the overall three-body Br is not important, so that 2b-3b interference should 
have little impact on the relative t]' vs. rj yields. However, hierarchies of 2 to 4 between 
the 7]' and rj yields may result from reasonable differences between the B ^ rj and B rj' 



^^It is also clear that it is not possible to obtain an 0(10%) charmless Br via new physics 
contributions to the penguin operators without violating rare B decay constraints, since this would 
necessarily increase the differential i?r's in the factorization region by an order of magnitude or 
more relative to the SM. For example, B{B — > (jjXg) would exceed the CLEO upper limit by a 
factor of 20 - 40 in the factorization approach. 
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form factors. In the case of enhanced b — > sg, Fig. (j^), the relative contribution of this 
region to the three-body yield appears to be significant, particularly for t] production, due 
to the singular behaviour at low g^, so that the interference could turn out to be relevant. 
For now only the pure two-body and three-body factorization estimates will be considered, 
with a low cutoff imposed on the latter. We make no pretense that factorization is valid 
at low for dipole operator contributions, recognizing that this further contributes to the 
uncertainty in t]' production in models with large b sg. 

In Table || we study the dependence of B{B^ rj'su) with experimental cut on choice 
of low cutoff. The singular behavior is apparent below 1 GeV^ for enhanced b ^ sg, 
while there is relatively little change in the SM Br^s below 1 GeV^. The three-body B ^ rj 
Sr's are very close to the rj' yields in the SM and ~ 5% — 15% larger with enhanced b — sg. 
In the ensuing discussion we fix the cutoff at 1 GeV^ for both t]' and r] production. We 
should note that the three-body r]^''^ yields are not very sensitive to m^, unlike the two-body 
yields. 

In Table Q the two-body and three-body CP averaged SM Br's are summed for the 
range of parameters previously discussed. We find that B{B v'-^s) lies in the range 
(0.5 — 2.5) X 10 which is significantly smaller than the observed rate even for the largest 
5r's. The CP asymmetries range from 2% to 5% with the larger asymmetries corresponding 
to the smaller i?r's. In the case of rj production the SM Br sums are slightly larger, lying 
in the range B{B r]X,) ^ (0.7 - 2.6) x 10-^ with CP asymmetries of 1.5% to 4.5%. 
As already noted, we do not expect significant changes in the ratio of r] to r]' yields when 
SM 2b-3b interference is taken into account. Finally, with no momentum cut we obtain 
B{B r]'Xs) ^ (0.25 - 1.0) x lO^^ in the SM. 

Taking B{b —>■ sg) ~ 15%, the dependence of rj' yields on the relative magnitudes of 
Cii and c'li, corresponding phases, and factorization model parameters is explored in Fig. 
5. Sums of pure two-body and pure three-body contributions have been plotted for three 
illustrative cases: = 0, Cn = 0, and |cii| = \c[-^\. It is worth noting that in models of 
quark mass generation a large hierarchy between cn and c'n is not expected in the absence of 
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some special flavor symmetry. According to Fig. (^) factorization alone could well account 
for the observed inclusive yield in Eq. (Q) if dipoles were enhanced.0 This conclusion persists 
when constraints from other rare decays are taken into account, as we'll see shortly. The 
large range of v( yields possible is in marked contrast to the low SM range, shown in grey. 
The Vj Br with no cut is ~ (2.5 — 10) x 10"^ for maximal constructive interference between 
dipole and penguin contributions. The r] yields are about 10% larger for F^''_(0) = -F-[L(0) 
and 2b-3b interference not taken into account. 

The CP asymmetries for inclusive 77' production in Fig. (|^) are not much larger than the 
SM range from Table |T|, despite the presence of significant new amplitudes with arbitrary 
weak phases. This is because the dominant three-body dipole contributions are weighted 
towards low g^, where little or no strong phase is generated for the interfering three-body 
penguin amplitudes. 



B. Exclusive ij-'^ production. 

It is straightforward to obtain factorization amplitudes for if-'^K^ from -f^e// and 

Hi\ in terms of the r^-'^ hadronic matrix elements listed in Eqs. (|13|, |l^, ^), the kaon matrix 
elements 

(ir-|s7^75M|0) = -zv^F^< (19) 

2 

(K-|s75m|0) = -i^PiFk ""I^" , (20) 



and the analogs of the B v( matrix elements in Eq. (|T^ for i? K transitions 
35|. Again we consider the range -Fi-(O) = .25 — .5 for the B — > r/*^'^ and B ^ K form 



factors at zero momentum transfer. Although we are looking at exclusive decays, individual 
quark momenta enter the amplitudes because of Hn. For simphcity we will identify ph with 



iipor F^_{0) = .5 but = .15 GeV, /§ = O^jn' « —17° the maximum yields would be about 
15% lower than in Fig. (^). 
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the B meson momentum and make the usual quark model kinematical assumptions, taking 
Ps = Ps = Prfof^ for non-spectator r/W production and Ps = Pu = Pk^/'^ for spectator rj^'^ 
production. With these assumptions one obtains = m^_/2 — m^(,)/4 + m|^_/2 for the 
former and = m^_/2 — m|._/4 + m^(,)/2 for the latter. 

SM Br's and CP asymmetries for rj'K^ are summarized in Table The Br's 

range from (1.1 — 12.1) x 10^^, with significant enhancement occurring at low values of 
due to the l/rris factors in (?7'|s75s|0) and {K~\s''y5u\0) . It is clear that SM factorization can 
reproduce the measured Br in Eq. (|^) given a reasonably broad range of input parameters. 
The CP asymmetries can also be significant, varying from 4.5% - 8 %. For rj production 
we find B{B^ r]K^) ^ (1.5 - 4.8) x 10"'' for F^'_{0) = F^_{0), in accord with the order 
of magnitude suppression expected relative to the rj'K^ yield from Lipkin's interference 
mechanism ||2^. The corresponding CP asymmetries are large, at the 15% - 20% level. 

Enhanced dipole operators can lead to significant increases in the exclusive rj' yield, as is 
to be expected given the increases we saw are possible for the semi-inclusive yield. Eq. (^ 
can therefore be reproduced for a larger region of factorization model parameter space. This 
is illustrated in Fig. (^), where exclusive t]' production with enhanced dipoles included is 
compared to SM yields for intermediate choices of form factors and m^. According to Fig. (^ 
significantly larger CP asymmetries than in the SM are also possible, although differences 
in magnitude decrease in regions of On and/or corresponding to large rj' production. 
Finally, it is important to note that Lipkin's interference mechanism for exclusive decays 
persists in the case of enhanced b ^ sg. For example, the rj yields corresponding to Fig. (^ 
lie in the range B{B r]K^) ~ (1.0 - 8.0) x 10"^ if F^'^{0) = F^^{0). 



C. Comparison with other rare B decays. 

It is important to check that large factorization model rj' yields obtained with enhanced 
b ^ sg are consistent with constraints from other rare decays. We have considered three 
processes: B (pX^, which is the most restrictive semi-inclusive decay, B^ — » (pK^, which 
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is probably the most restrictive exclusive decay, and K'^ir^ which is not particularly 

restrictive but, like the other two, is least cumbersome to calculate because there are no tree- 
level contributions (pure-penguin decays). Recently the CLEO collaboration has presented 
the following bounds: B{B (pX^) < 2.2x10-^ (90% c.l.) |2|], B{B^ <pK^) < 1.2x10"^ 
(90% c.l.) 0], and B{B^ K^tt^) = 2.31:\m x IQ-^ In the SM we obtain B{B ^ 
(pXs) ^ 1.1 X 10-^ (also see @ii]), and B{B^ ^K^) ^ (0.32 - 1.3) x 10"^ as F/^(0) 
is varied from .25 to .5. Neither Br is sensitive to m^. For F[-{0) = .333 |^ we obtain 



B{B^ — > K^Ti^) = (1.0 — 2.1) X 10 ^ as is varied from .2 to .1 GeV. It is interesting that 
the largest SM rj' yields in Tables |T| and |V| are consistent with the above constraints. For 



the exclusive decays this is due to Lipkin's mechanism, whereas in the semi-inclusive case this 
is because there is no spectator contribution to production, whereas it constitutes the bulk 
of the f]' yield. Another semi-inclusive process for which there is a spectator contribution, 
B KsX, is very weakly bounded, B{B K,X) < 7.5 x 10^^ for pk > 2 GeV [gg. 

In Fig. d^) we show the (pXs and (pK^ yields with enhanced b ^ sg corresponding to the 
examples in Fig. (^. Comparison of Figs. (|§ - 0) shows that within the factorization model 
there is significant overlap of regions in (cn, c'n) space with exclusive and semi-inclusive r]' 
yields in the measured ranges and acceptable (p production. No additional constraints arise 
from B^ K^TT^. 

III. EXOTIC MECHANISMS 

It is well known that the 1]' is not a Goldstone boson - its mass does not vanish in the 
chiral limit mu^d,s 0. This is believed to be due to a large glue component in the rj'. In 
the present context the enhanced glue content leads to a significant coupling of the rj' to two 
gluons, which participates in 6 — > s transitions. 

It is important to notice that because of large cancellations p!9| , ^ , the contribution of 



the b sgg mode followed hj a, gg ^ rj' transition is extremely small. On the other hand, 
the (9(1%) Br of the b sg* transition in the SM could, in principle, be responsible for the 
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anomalously large b rj'Xg decay rate. This fact was exploited by AS [T^, who considered 
the contribution from the h sg* rj'sg subprocess. They parameterized the effective 
ggrj' coupling as 

V;,e^^ = H{q\ k\ ql,)e^p,.q^k^e^,e'^ (21) 

where = {pb —PsY- A constant form factor was assumed, i.e., H {q"^ , k"^ , q"^,) ~ 
H[0,0,m^,) = Hq, and Hq was extracted directly from the decay rate for J/ip ^ T]''y. 
For as{mj/^) = 0.25 and Orjrt' = —17°, the value Hq = 1.8 GeV^^ is obtained. 

With the central assumption of weak dependence in H, AS found that the observed 
7]' yield could be fully accounted for. However, it is clear that in order to obtain the total 
decay rate, the differential distribution for this subprocess must be integrated over a wide 
range of g^, spanning approximately 1 GeV^ to m^. It is therefore of paramount importance 
to investigate the off-shell dependence of the form-factor describing the g* — * rj'g transition. 

Determination of the exact g^-dependence of H{q^, 0, Tn^, ) is a difficult task. We therefore 
employ a semi-phenomenological description of the form-factor, picking up the leading g^- 
dependence and parametrizing the rest by an appropriate constant. To set the stage we 
recall the leading g^-dependence of the 7* Vj^ form-factor in the limit of large momentum 



transfer ^l|. In this limit, a perturbative QCD description is possible and it is not difficult 
to convince oneself that the simplest two-particle irreducible kernel yields a form-factor that 
scales like l/Q^ as ^ cx) (Q^ = _g2) |||. 



r^^e^^ = ze2F,,(Q2)e„/3^,g'^A;^e^^, F^,{(f) = 4^^^^- (22) 

with ^ ~ 0.5. While the corresponding calculation for g* — ^ r]'g is more involved, the lowest 
order dependence in the — > 00 limit should be the same, so it is clear that the constant 
form factor approximation H{q^, 0, m,^/) ~ Hq is not acceptable. 

In order to obtain a more realistic parameterization we consider a model of the ggr)' 
vertex in which we couple a pseudoscalar current to two gluons through quark loops. This 
approximation is believed to be a good one since the "direct" ggrj' coupling is suppressed - 
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the quantum numbers of t]' are , while the lowest possible bound state of two gluons is 
0"*". The calculation yields 



V, 



E 



afgtmf 



rl rl—x 

/=M,s 27r2 Jo Jo 



dy 



y- — (1 — X — y){xq^ + yk"^) — xym^i — ie 



m 



E 

f=u,d,s 



f 



^-^F^I\q')e.,,.q-k^ F'i\q') = I {^) - I . (23) 
-mi, ^mi^ ^mV 



27r2 g2 



where the /(a;^)-functions are defined as 



-2arcsin 



M 

2 ' 



< < 4, 



(24) 



a/ is the coupling of the t]' to the quarks in the loop, i.e., afrj'f'y^f. One has to keep in mind 
that this vertex is in general non-local, so that aj will depend on the momentum transferred 
through the vertex. Since this dependence is at most logarithmic, i.e. subleading, it will be 
ignored in what follows. 

Motivated by the above result it is now reasonable to parametrize the effective ggrj' 
coupling as 



V, 



Homl, 



\(^ai3f,uq°'k'^, and H{q^,0,ml,) 



Hom^, 
— mil 



(25) 



Note that Ho does depend on q^. We drop this dependence based on the following observa- 
tions. First of all, q^ > Am'j for any quark flavor over the whole range of contributing to 
the total decay rate. This implies that F^\q'^) depends on only logarithmically. Also, as 
noted by HT |16[, the factor Q;s(/i^) implicit in Hq must in general be running. They argued 
that the scale at which is evaluated should be associated with the momentum transfered 
through the ggr]' vertex, introducing another logarithmic dependence on in Ho- It is clear 
that these dependences are subleading with respect to the strong power dependence in Eq. 
(p5|) and to first approximation can be modeled by a constant Hq, which we identify with 
the value extracted from J/ip —>■ r]''-/. 

The strong g^ dependence of Eq. ( P^D implies additional suppression of the b rj'sg 
anomaly subprocess at large g^. It is easy to estimate the resulting contribution to 6 ^ rj'Xs- 
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The calculation is a straightforward generalization of the one carried out by AS [T^ and HT 
p!6| so we do not include any details here. We only point out that in the SM the interference 
between the charge radius and dipole form factor contributions is destructive. 

Our SM result for the total branching ratio is B{h —>■ rj'sg) ^ 5.3 x 10~^, corresponding 
to B{b — s> rj'sg) ^ 1.6 x 10^^ with experimental cut. The latter is more than an order of 
magnitude smaller than the observed t]' yield in Eq. The recoil spectrum is shown in Fig. 
8. As already pointed out in Refs. |lT5| , p!6| , the anomaly subprocess exhibits a rise with 
in the low mx^ region which is characteristic of three-body decays. The dependence of 
the ggrj' formfactor significantly distorts the distribution so that the maximum is shifted 
towards lower values of compared to constant ggt]' formfactor. 

In the case of enhanced b ^ sg we again consider B{b sg) ~ 15%, and for simplicity 
take c'li = 0. The resulting Br^s with experimental cut are presented in Fig. 9 as a funtion 
of 6[i. These are more than an order of magnitude smaller than obtained by HT under 
the assumption of mild g^-dependence of the form factor. The corresponding CP violating 
asymmetries are also systematically smaller than those obtained by HT [|T^ . This is because 
of the shift in the decay distribution towards lower values of g^, leading to smaller strong 
phases in the charge radius form factor. 

We also remark that the short- distance b rj'sg subprocess most probably does not 
affect the exclusive B —>■ rj'K branching ratios. If the gluon attaches itself to the s quark 
then the exclusive decay can be interpreted as being due to the short distance subprocess 
b — i> sgg which we already know is suppressed [|19|. The other two possibilities, that the 
gluon attaches itself to the spectator quark, or that it contributes to the process via higher 
Fock states of the participating mesons, must be additionally suppressed by the dynamics 
of the process. 

Another interesting explanation of the anomalously large r]' production rate is to assume 
a large 'intrinsic charm' content for the rj' Phenomenologically, this can be described 
as a mixing of the rjc and rj' . Since the value of the mixing angle is obtained by a fit to the 
experimental branching ratios for radiative charmonia decays to rj' this phenomenological 



estimate should already contain a contribution from the higher Fock states containing charm 
quarks [jl8[. It is worth mentioning that intrinsic heavy quarks of the light mesons must 
be far from their mass shells in order to satisfy energy-momentum conservation. Therefore, 
NRQCD methods developed for weakly bound quarks in charmonia can not be used to 
describe the dynamics of intrinsic heavy quark states. 



The heavy-light quark mixing for the rj—rj'—rjc system was considered in [^] and found to 
be small, sinap ~ 2.4 x 10"^. This leads to the estimate Br{B v'-^s) ~ sin^ apBr{B — >■ 
T]cXs) < 10^^ so that the 'intrinsic charm' content of the t]' is not likely to be responsible for 
the observed large branching ratio. In addition to phenomenological estimates limiting the 
size of the cc contribution, a large 'intrinsic charm' explanation for the rj' yields would lead 
to B{B — s> rj'K*) ~ 2B{B — > r]'K) and a quasi two-body momentum spectrum for B VjXg 
[T7|, neither of which is observed. 



IV. CONCLUSIONS 

We have considered production of v( in B meson decays in the Standard Model and in 
models with enhanced h — sg. It is clear that SM factorization can, in principle, account 
for the exclusive Vj yield. For a liberal range of factorization model parameters the SM 
B^ rj'K^ branching ratio lies in the range (1.1 — 5.8) x 10~^ for = .2 GeV and 
(2.3 — 12.1) X 10"^ for = .1 GeV . Sizable CP violating asymmetries at the 5% - 8% level 
are expected. 

The inclusive yield requires a more complicated solution. For the same range of param- 
eters SM factorization gives B{B'^ v'-^s) ~ (0.5 — 2.5) x 10~^ including the experimental 
cut, with the largest yields corresponding to a fairly limited region of parameter space. In 
addition, the branching ratio of the QCD anomaly mediated subprocess b rj'sg considered 
in [|I5| is reduced to ~ 1.6 x 10~^ when taking proper account of the leading g^-dependence 
of the rj'gg coupling form factor, H{q^) ~ m^i / {q^ —mf^,). The large rate reported in |]T5| for 
constant form factor would require a factor of 50 enhancement due to non-perturbative ef- 
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fects, e.g., resonances, of the rate obtained with leading g^-dependence included - an unlikely 
possibility. However, a SM 'cocktail' solution for large B{B^ v'-^s) involving contribu- 
tions from several mechanisms, e.g., factorization pushed to the limit (the largest single 
contribution), large non-factorizable contributions, charmonia decays [|1^], 'intrinsic charm' 
]r7| , and b rj'sg, is still a possibility. 

On the other hand, the intervention of Non-SM physics in the form of enhanced chromo- 
magnetic dipole operators provides a simple and elegant solution to the puzzle of large rj' 
production in B decays. This explanation can be realized in a broad region of the available 
parameter space. In fact, factorization model contributions alone can account for both the 
semi-inclusive and exclusive rj' yields without violating contraints from other rare decays. 
The b rj'sg rate, although significantly smaller, can be an order of magnitude larger than 
in the SM. CP violating asymmetries in exclusive and semi- inclusive B ^ rj' decays due 
to new weak phases in the dipole operator coefficients can be larger than in the SM and 
of either sign. However, large asymmetries in the 'pure-penguin' decays B"^ (pK^, K^n^ 
would provide a more definitive signal for new physics because of the small asymmetries 
expected for these decays in the SM. 
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TABLES 

TABLE I. Two-body contributions to B{B^ — ^ v'^s) (average of CP conjugate decays) and 
corresponding CP asymmetries for (i) standard model without On, (ii) standard model with On. 
-^0 = Ft^- 



Inputs 


no On 


with On 


ms = .1, = F,, Ortn' = -17° 
ms = .2, Fs = F^, Orjrj' = -17° 
ms = .1, Fs = 1.25F^, 0^^/ = -21° 
ms = .2, Fs = 1.25F^, 0^^, = -21° 


3.7 X 10-^ .51% 

1.3 X 10-^ 1.5% 

8.4 X 10-^ .21% 

2.8 X 10-^ .63% 


3.2 X 10-^ .59% 
1.04 X 10-^ 1.8% 
7.4 X 10-^ .24% 

2.34 X 10-^ .75% 



TABLE IL Three-body contributions to B{B v'^s), with experimental cut, versus 
[GeV'^] for (i) standard model with On, (ii) B{b sg) 15% with |cn| = |c'n|, 9n = 180°, 
^ii = 0. F;''(0) = .5, = .1 GeV. 



o2 . 


SM 


B{b sg) ^ 15% 


.4 


1.71 X 10-4 


12.1 X 10-4 


.6 


1.68 X lO"'' 


9.61 X 10-4 


.8 


1.66 X 10-4 


8.38 X 10-4 


1.0 


1.65 X 10-4 


7.61 X 10-4 


1.4 


1.64 X 10-4 


6.64 X 10-4 


2.0 


1.62 X 10-4 


5.77 X 10-4 
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TABLE III. Sum of two-body and three-body contributions to B{B^ v'^s) (average of CP 
conjugate decays) in SM for (i) F^'_{0) = .25, (ii) F^'_{0) = .5. Fo = F^. 



Inputs 


F^'_{0) = 


.25 


F^''!(0) = .50 


"ms = -1, -^8 = Ft^, dr)-n' = 


-17° 


7.6 X 10" 


-5 


2.1 X IQ--^ 


f^s = -2, Fg = Ftt, Or^ri' = 


-17° 


5.4 X 10- 


-5 


1.8 X 10-4 


ms = .1, Fs = 1.25F7r, 6*^^' 


= -21° 


1.2 X 10" 


~4 


2.5 X 10""^ 


nis = .2, Fg = 1.25F^, 


= -21° 


6.7 X 10" 


-5 


2.0 X 10-4 



TABLE IV. B{B^ rj' K^) ( average of CP conjugate decays) and corresponding CP asym- 
metries [%] in SM. Fo = F^. 



Inputs 


rus = .1 GeV 


nis = -2 GeV 


= F^, Or^rj' = -17°, Fi'''(0) = Ff (0) = .25 
Fs = F„, e^^, = -17°, Fi^'lO) = Ff (0) = .5 
Fs = 1.25F^, er„: = -21°, Fi'''(0) = Ff (0) = .25 
Fs = 1.25F^, e^r,' = -21°, Fi'''(0) = Ff (0) = .5 


2.3 X 10-^ 5.4% 

9.1 X 10-^ 5.4% 
3.0 X 10-^ 4.5% 

1.2 X 10-4, 4.5% 


1.1 X 10-^ 8.0% 
4.6 X 10-^ 8.0% 
1.5 X 10-^ 6.9% 
5.8 X 10-^ 6.9 % 
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FIGURES 
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FIG. 1. W^^{b sdd) for: (a) SM, (b) Enhanced On by itself, (c) SM with enhanced On 
for Oil = 180° (constructive interference). The regions relevant to two-body or quasi two-body 
decays in the factorization approximation are blown up. 
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FIG. 4. 10^3^(6 ^ sdd) for: (a) SM, (b) SM with enhanced On and = 180°, in the 
regions relevant to ij' production with experimental cut. 
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FIG. 5. B{B^ rj'su) vs. 9 and CP asymmetries vs. 9 for B{b sg) ~ 15%. Shaded regions 
are the corresponding SM ranges. Sohd curves are for c'n = {9 = 9ii), long dashed curves are 
for |cii| = \c'ii\ and 9'ii = 180" {9 = 9ii), short dashed curves are for cn = {9 = 9'ii). The larger 
Sr's (lower asymmetries) are for F^_(0) = .5, /g = 1.25/,^, 9^^^' = —21°, m<j = .1 GeV, the lower 



Sr's (larger asymmetries) are for F^_{0) = .25, /§ = f^, ^r^rj' = —17°, m<j = .2 Gel^. 
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FIG. 6. (a) 10^B{B^ rj'K^) vs. 6 and SM value (horizontal line) for F'^_{{)) = .35, 
F^_{Q) = .38, ms = .15 GeV, Fg = F^, Onn' = -17°. Solid, long dashed, and short dashed 
curves are for same cn and c/n as in Fig. 5. (b) CP asymmetries vs. 9 corresponding to (a). 
Horizontal line is SM value, (c) B{B^ — r]'K^) in (^11,^11) plane for |cii| = \c'ii\. 
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FIG. 7. (a) W^B{B ^X^) vs. 6 for 13{b sg) w 15%. Solid curve is for c'n = {9 = 0u), 
dashed curve is for cn = {9 = 9[i), horizontal solid curve is SM value, and dot-dashed curve is 
the CLEO upper limit, (b) B{B — > (pXg) in the {9ii,9'ii) plane for |cii| = \c'ii\. (c) and (d) are 
analogous to (a) and (b) for 10^^(5=^ (t>K^)- In (c) solid and dashed curves coincide. 
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FIG. 8. ^^-(6 —>■ vs. in the SM. Curves from top to bottom are for rrn, = 5.2, 5.0, 
and 4.8 GeV. nig = 0.5 GeV and = 0.2 GeV. 
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